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Abstract 

We consider the production of open charm in heavy ion cohisions in the frame- 
work of the Color Glass Condensate. In the central rapidity region at RHIC, for the 
charm quark yield we expect A'^cou (number of collisions) scaling in the absence of 
final-state effects. At higher energies, or forward rapidities at RHIC, the saturation 
scale exceeds the charm quark mass; we find that this results in the approximate 
A^part (number of participants) scaling of charm production in AA collisions and 
\/^part scaling in p{d)A collisions, similarly to the production of high pT gluons 
discussed earlier. We also show that the saturation phenomenon makes spectra 
harder as compared to the naive parton model approach. We then discuss the en- 
ergy loss of charm quarks in hot and cold media and argue that the hardness of the 
spectrum implies very slow dependence of the quenching factor on pT- 

1 Introduction 

Heavy quarks play a very important role in the development of Quantum Chromody- 
namics. The masses of heavy quarks m significantly exceed the QCD scale Aqcd, which 
makes perturbative calculations of charm production and annihilation possible PJ. The 
QCD decoupling theorems |2j ensure that heavy quarks do not influence the dynamics of 
processes at scales much smaller than m. 

At sufficiently small Bjorken x and/or for sufficiently heavy nucleus A, parton distri- 
butions approach saturation 0. In the saturation regime partons form the Color Glass 



Condensate (CGC) characterized by the dimensionful "saturation scale" Qs determined 
by the parton density [SI IH El E] ■ Possible manifestations of the Color Glass Condensate 
in the energy, centrality, and rapidity dependence of hadron multiplicities [3 |Hl E] have 
been found by experiments at RHIC [IHl HD IH CHI • 

It is intuitively clear that once the parton density becomes high enough to ensure that 
Qs and m are of the same order, heavy quarks will no longer decouple. This leads to very 
interesting consequences for their production dynamics, as we will discuss at some length 
in this paper. First, let us however illustrate this point by using very simple qualitative 
arguments. 

The Color Glass Condensate is characterized by strong classical color fields [HI El ; the 
strength of the chromo-electric field can be estimated as 

E-^, (1) 
9 

where = Airag is the strong coupling constant. This field can polarize the Dirac vacuum 
of quarks; in particular, quark pairs can be produced when the potential difference across 
the Compton wavelength of the heavy quark becomes equal to the energy needed for the 
pair production, 

m r, 

gEr^-— = m\ 2 

1/m 

Together with (P) this condition suggests that when 

Ql > m\ (3) 

the heavy quarks will no longer decouple, and their production pattern will be similar to 
that of light quarks. 

It has been argued before that Color Glass Condensate at sufficiently high energy 
and/or rapidity leads to the suppression of high pt parton production. Namely, instead 
of the Ncoii (number of collisions) scaling of yields of high pt partons expected on the 
basis of QCD collinear factorization, one finds Npart (number of participants) scaling in 



AA collisions and y A^^j^rt scaling in p{d)A collisions. In view of the argument given above, 
one expects that the yields of heavy quarks can follow similar scaling. The goal of the 
present paper is to establish in what kinematical region this scaling can hold for heavy 
quark production. 

The present d-Au data [221 mid-rapidity at RHIC do not support the suggestion ^ 
that quantum evolution in the CGC is responsible for the observed in Au — Au collisions 
suppression of high px hadrons. Instead, they indicate the dominance of final state effects 
consistent with jet quenching [121 1201 1211 • Nevertheless, at sufficiently large energy and/or 
rapidity the arguments of [HI should apply. Since the saturation scale exponentially in- 
creases with rapidity, one may expect that this can happen already at RHIC energies in 
the forward rapidity region. 

Production of particles in relativistic heavy-ion collision is a complicated process which 
in the CGC framework can be considered as consisting of two stages. First, upon the 
collision partons are released from the wave functions of the nuclei. At the quasi-classical 
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level the wave function of the ultra-relativistic nucleus is described by the non-abelian 
Weizsacker- Williams field [HI ^5 created by the valence quarks moving along the light cone. 
Quantum evolution of the wave function of the nucleus is described in the framework of the 
Color Glass Condensate by a set of coupled evolution equations [12] which can be reduced 
to a single nonlinear equation in large approximation ^H]- At the second stage, there 
are final state interactions between the produced partons which may bring the quark- 
gluon system to thermal equilibrium PJj- In this paper we address both problems in the 
case of the open charm production. Very recently, heavy quark production in the CGC 
framework has been considered also in Ref . jTH| . Our treatment of the problem includes the 
effects of quantum evolution in the nuclear wave functions; they appear to be significant 
and lead to non-trivial centrality dependence of the heavy quark yields. 

2 Heavy quark production in /c^-factorization 
2.1 Collinear and kx factorization schemes 

Unlike the inclusive gluon production case 0123 El El, as we discussed above the heavy 
quark production is characterized by two inherent scales: the saturation scale Qs and 
the heavy quark mass m (we assume that the two colliding nuclei have the same atomic 
number A). Depending on the relation between these two scales we have two different 
mechanisms of the quark production. The quasi-classical gluon field of a nucleus is char- 
acterized be the unintegrated gluon distribution fA{k±,y) which is fiat at k'j_ < Ql* and 
decreases as l/k\ otherwise. A quark production amplitude A in the process gg qqX 
has typical momentum p\ ~ max{m^,Q^}. The cross section of the quark production 
is given by the convolution of the unintegrated gluon distributions of nuclei with the 
quark production amplitude. This statement is known as kr factorization [231l2nil2Zl- It 
was proved only in the case when multiple rescatterings of partons from the nuclei wave 
functions are suppressed relative to the evolution effects. Although the fcy factorization 
has not been proven for general process in high energy QCD, its numerous applications 
are proved to be rather successful (for charm production, see e.g. j2Hl)- In the context 
of relativistic heavy ion program it was utilized in jH[ E] to calculate the inclusive gluon 
production. Additional confidence in /cT-factorization arises from the fact that it holds 
exactly (includes all rescatterings in a nucleus) for the gluon production in 7*/! [29] . 
This result also holds for a pA process at not very high energies, when a proton can be 
considered as a diluted object [21]. In the framework of the CGC fc^ factorization for 
the heavy quark production was recently demonstrated in ^H]- In the following we will 
assume that the fcr-factorization gives fairly good approximation at RHIC energies in the 
central rapidity region even if the rescatterings are turned in. 

The basic assumption of the collinear factorization is that the typical transverse mo- 
mentum associated with the colliding hadron is of the order of the Aqcd and much smaller 
than the typical momentum of the hard subprocess. In the case of heavy flavor production, 
large quark mass m seems to insure the validity of the collinear factorization. However, 
at high energy the typical scale of the hadron wave function is Qs{y) which increases 

*See extensive discussion in |24| . 
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with energy. It becomes no longer possible to integrate out the transverse degrees of 
freedom of the hadron wave function. The /cT-factorization is thus the generalization of 
the collinear factorization at high energies. When Qg '^m a quark production amplitude 
overlaps only with perturbative {pt ^ Qs) low density tail of ip. Perturbative approach 
must be valid in that kinematic region although we will comment later how strong must 
the above inequality be. In the opposite case Qs ^ m unintegrated gluon distributions 
and the quark production amplitude strongly overlap. The hard transverse momentum of 
gluons in a nucleus field can be transmitted to the quark-anti-quark system giving rise to 
a strong deviation from the naive perturbative approach. In particular, the high parton 
density effects might manifest themselves in the form of the charmed-meson spectrum, 
dependence of the quark multiplicity on the atomic weight A, angular asymmetry of quark 
jets production and overall enhancement of the quark production cross section. 



2.2 A^part versus A^coii scaling in inclusive gluon production 

The question about the relation between scales Qs and m is closely related to the question 
about the dependence of the quark multiplicity on the atomic weight A. Let us first review 
the argument given for the inclusive hadron production case jHl El 121] • The multiplicity 
spectrum is given by 

J dklifAixi, kl) <fA{x2, (P - k)']_). (4) 



dN 27ia, 1 



dyd'^p±_ Cf Sa p]_ 

We would like to consider it in three kinematic regions which are defined by three di- 
mensional parameters Qo? Qs and fcgeom = Q^/Qo- Here Qq is the characteristic scale 
of nuclear hadronization region. The scale /Cgeom sets the region of applicability of the 
collinear factorization. Indeed, at k± <^ fcgeom we have (modulo logarithms) 

Inkl/Ql = ^^kl/Ql + \nQl/Ql ^ \nQl/Ql = Xy, (5) 

where A ~ 0{as)- Thus, the evolution with respect to the longitudinal momentum 
(BFKL) cannot be neglected with respect to the evolution in the transverse momen- 
tum (DGLAP) despite the fact that evolution is linear. We will estimate the integral 
in ^ assuming that k'j_ -C p^. This corresponds to the logarithmic (in p±) approxima- 
tion. Note, that integration in (@]) goes always over the smallest transverse momentum 
min{k\, {p — k)^}, while the external one is the largest. This can be easily seen by rewrit- 
ing in symmetric from in terms of k_i_ and qj_ = p± — kj_. In the region of applicability 
of the parton model p\ > fcgeom we have 

dN Sa fPl,,2Q'sQ's Q n4 1 at ff,^ 
dkj^jY^ SaQs T ~ ^cou, (6) 



dyd'^pi_ asp\ J k\ p\ " as pi 

smce 

Q2 ^ ^1/3 and Sa ~ A^/\ In the region Ql < p\ ^geom evolution is linear but 
existence of the strong color field at long distances affects the gluon distribution which 
acquires the anomalous dimension 1/2: f ^ k±/Qs |SD1ISI1E2]- Thus, 



which is known as the geometric scahng . Finally, in the saturation region p]_ <^ Ql 



dN Sa 



dki ~ Saq: — r ~ ^p-t (8) 



We reproduce the result of [H] that the multiplicity scale with A^part in a vast kinematic 



region k± < fcgeom- Similar considerations lead to y Np^rt scaling in p(d)A collisions for 
sufficiently central collisions. 

The transition from the classical regime dominated by the Cronin effect |^ to the 
suppression regime was discussed in our recent paper with Yu. Kovchegov [25. This 
issue was addressed recently in Refs. [HI], which confirmed the phenomenon of high px 
suppression at sufficiently large rapidity and studied its evolution. 



2.3 Heavy quark production 

Returning to the case of heavy quark production, the saturation corresponds to the region 
^± < Qs, where = p^ + rn^ ■ To estimate the scale of the coUinear factorization 
breakdown we note that the heavy quark threshold in the s-channel is 4m^. Analogously 
to the gluon production case considered above, the scale at which collinear factorization 
breaks down is 

M = % (9) 

where M. is the invariant mass of the heavy quark pair; for the quarks produced at close 
rapidities, A^^ ~ 4 {j>\ + rn?) (in their CM system). Since at RHIC Qg is of the same 
order as m in the midrapidity, is not satisfied even at p^ = 0. Therefore, we expect to 
observe the Ncoii scaling of the open charm spectrum. At the forward rapidity Qg ^ m 
and we expect to observe A^part scaling for momenta p]_ < max{Q^, Q^/(4(5o) — rn?}- In 
any case, the charmed-meson spectrum is expected to be much harder than that predicted 
by the parton model, since the typical gluon momentum in the nucleus wave function is 
Qs 3> Aqcd- The total cross section is also expected to be higher since there are additional 
contributions to the production amplitude as compared to the parton model. 

It was argued in Refs. [23 123 12Z| that the cross section of the hadroproduction of the 
heavy quarks at high energies can be written in a /c^-factorized form: 

da f d'^q±i r d'^pj_2 



/ / ^A{y,q±i)VA{x,{p±i+p±2- q±i) ) 

J IT J 71 



Here p^, p^ and yl, are the quark and antiquark transverse momenta and rapidities, 

S = {pa +PBf and 

^1,2 = {<i,2lSf'^e-yl^, X = X,+X2 (11) 

2/1,2 = «i,2/5)i/2g,*,,^ y = yi + 2/2, (12) 
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Figure 1: Production of qq pair in heavy- ion collisions in fc^-factorization approach. 



q\) is unintegrated gluon distribution in a proton defined as 

2^ dxGA{x,qi) . 

fA{x,qJ = . 13 

dqi 

A{gg ggcc) is a production amplitude the exact expression for which is rather nasty 
and can be found in the Appendix. In the limit of vanishing virtualities of the t-channel 
gluons g^, ~^ equation (fTUj) reduces to the well known parton model expression 
(collinear factorization) 

The amplitude ()47j) is rather hard for analytical analysis. We can simplify it in the 
Regge kinematics ly* — I/2I ^ 1- This limit cannot be used in numerical calculations since 
quarks tend to be produced at the same rapidity: production of large invariant masses 
of the qq pair is suppressed. However it gives a fair qualitative understanding of the 
underlying physics. A simple calculation yields 

. A{gg ^ ggcc) ^ + 0{e-'^y ). (14) 



In the region of applicability of collinear factorization, mix 00, which at RHIC is the 
same as mix ^ Qsmax, the transverse momenta of gluons are negligible compared to those 
of produced heavy quarks. Therefore pix ~ —P21. and (fTUj) can be written as (we omit 
numerical factors) 

dN""^ ale-^y* fPl± dpi, /■ , 2 1 A^xGN(qh)xGN(pli) 
(X — / / dgi 



d^p^idyldyl Sa J mi^max 7 A^cP?_l 9i±P2± 

alA^e-^^y^xG^ipl^))^ 



NcSap\^ 



A^coii, (15) 



where clearly p\^ ^ m?. Note that the exponential suppression at large Ay* of the heavy 
quark amplitude as compared to the gg — > gggg [SEl is because t-channel quark carries 
spin 1/2. 
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At large rapidities there is a kinematic region such that Qsmm ^ '^li ^ Qsmax- In 
this case one of the nuclear wave functions is in the saturation whereas another one is 
not. In general transverse momenta of quarks are no longer equal and point back-to-back. 
Thus in the limit pf^ ^ the quark multiplicity reads 

dN^' oc ^'e-^^* f dpl^ rQl^.. 1 Sa QLin Sa 

d^p±idyldy2 Sa J mi±m2± J ^"^^cP?_l 

^ iv p'^'"* ^ ^ 

^^c Fl_L Vsmin Vsmin 

Please, note that here Ql^^^ ~ asA/S^. 

Finally, there is a region where both nuclei are in saturation ^ Qsmin- We have 

dN^^ ^ ale-^y* r dpl^ rQL.n 1 Si 



d'^p±idyldy^ Sa J mi^max 7 N^plj^ a 

1 e~^y*SAQ^ 



/ smin 



iVparf (17) 



Consider now the total multiplicity of heavy quark pairs produced in heavy ion colli- 
sions per unit rapidity, dN'^'^ /dyldy2- At RHIC in the kinematic region |y*|<l we have 
= 0) ~ 1 — 2 GeV^. Since m ^ 1.2 — 1.9 GeV [HT] the transverse size of charmed 
quark 1/m is too small to feel saturation in either one of the nuclei and set A^part scaling. 
The parametric dependence of the total multiplicity of heavy quark pairs is (see ()15|) ) 

° -N^oii, 2/*^0, RHIC. (18) 



dyfdy^ S, 

Nevertheless, the fc^-factorization result for total multiplicity is numerically different from 
the parton model one: since in fc^-factorization spectrum is harder we get larger total cross 
section. 

Note, that since at RHIC Qs ^ m the variation of the saturation scale with rapidity 
becomes essential. Variation of Qs with rapidity stems from the fact that the saturation 
scale is proportional to the gluon density in the transverse plane of nucleus, which in turn 
is proportional to the gluon structure function xG{x,kj_). Since xG{x,k\) is increasing 
function of rapidity y* = ln(l/a;) we conclude that Qg grows with y* defined with respect 
to the rapidity of the nucleus in a given reference frame. Therefore, relation between 
saturation scales of each nucleus and the quark mass are different at different rapidities 
y*. If we consider quarks produced at \y*\ ^ 1 (but still not in hadronization region), 
then one of the nuclei is in saturation while another one is not. In that case we estimate 
using (fTHj) 

~ asA ~ ATpart, \y*\ > 1, RHIC, LHC. (19) 

"1/2 

In p{d)A collisions, we expect that in the deuteron fragmentation region and for sufficiently 
central events, the scaling will be 

dN'''^ 



dyldyl ^ \y*\ » 1' LHC. (20) 
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In the next section we will quantify the kinematical region at which the transition to the 
new scaling will take place at RHIC. 

At LHC energies the situation is different. We expect that the saturation scale in the 
midrapidity region will be much larger than the charmed quark mass. As the result both 
nuclei will be in the saturation. By (fT7|) we get 

— ~a«A~7V,„. LHC. (21) 

whereas far from midrapidity one of the nuclei is not in saturation and we can use (|19|). 
Note that from (|lH) . (jl2p and (|36|) it follows that the ratio between the largest and the 
smallest saturation scales is 

^smax ^ gA(|y*| + |yJ|)_ ^^22) 
Q smin 

This ratio is much larger than unity when at least one of the quarks is produced at 
\y*\ ^ 1 although in most cases quarks are produced at ll/* — y2l~l since otherwise the 
amplitude is exponentially suppressed (see (HH)). 



3 Open charm production in AA and dA collisions 

3.1 A model for unintegrated gluon distribution 

To proceed further we have to specify the unintegrated gluon distribution of a nucleus 
g). In principle it is directly related to the forward elastic scattering amplitude, 
which in turn can be calculated from the nonlinear evolution equations . However, 
their exact analytical solution is not known. Although recently the progress in obtain- 
ing the numerical solution was reported by many authors [HH] we prefer to use simple 
parametrization which catches the most essential details of the solution. In this paper we 
employ a model for distribution function suggested in 9J. It matches the known analytical 
expressions in the asymptotic regions [30,^. It reads 




Ql) 



(23) 



The modified Bessel function Jq is a solution of the DGLAP and BFKL equations in the 
double logarithmic approximation. 



.M^(^..^±i-l. ,24) 

In derivation of (j^Hj) the phenomenological parametrization of anomalous dimension in 
the Mellin moment variable uj was used j4Q] . 

= - l) . (25) 
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The first term in the right-hand-side of is the leading order contribution to the anoma- 
lous dimension of the DGLAP and BFKL equations in the double-logarithmic approxi- 
mation. The second term is the phenomenological correction which imposes momentum 
conservation 7(1) = and describes DIS data well. We also imposed the correct large x 
behavior on ^p^ by multiplying it by a factor (1 — x)^. 

Let us consider the unintegrated gluon distribution in various kinematic regions. First, 
consider a perturbative regime in which the collinear factorization is valid q±_ ^ fcgeom(z/*), 
see dSj) and 0- Performing expansion in the exponent of the Bessel function in this regime, 
using its well-known asymptotic behavior Iq{z) up to a slowly varying logarithmic 
factors and expanding in ()24|) (i(r) Ql{y*)/q]_ we get 




Next, when the typical parton momentum in the heavy- ion wave function approaches the 
saturation scale q± -C fcgeom(2/) one can expand the Bessel function as jH]: 



In In ^ 1 - In ^ H iQliy*)/<liy^^ (27) 




exp 

which yields 

(I--)- (28) 

Note that the anomalous dimension is 7 = 1/2 + 0{as) as desired. Finally, in the 
saturation region g_L < Qsijj) the unintegrated gluon distribution reads 

^a(^' ^ 4^2— HQVql) (1 - ^f- (29) 

2(27r)^ as 

This expression coincides with the formula for the unintegrated gluon distribution in the 
saturation regime in the quasi-classical approximation (2^. Actually, it also holds beyond 
the quasi-classical approximation. Indeed, (Pa{.x, q^) can be expressed via the gluon dipole 
forward scattering amplitude Nc^x^r) as [211 Elj 

VA{x,q^) = ^'^Z^.i I drr Jo{kr)\ NG{x,r^), (30) 
- as[2TT)^ Jo 

where r is the dipole transverse size. As was argued in jHOj, the forward scattering 
amplitude in the momentum space at x ^ and q fixed is 

N{x,q) = HQJq). (31) 

In the configuration space we have 

N{x,r) = [^^ dkkJo{kr) ln(g,/fc) = 1 - MQsr). (32) 
Jo 
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Using the unitarity constraint Nq = 2N — N"^ and we derive 




(33) 



(34) 



which coincides with fl29|) . 

Dependence of the saturation scale on energy is known in the double-logarithmic 
approximation |3(H l39] . It was shown in jUj that the photon structure function in the 
DIS in the whole kinematic region x < 0.01 can be described by simple Glauber-like 
parametrization of the gluon distribution with the following saturation scale 



where xq and A are a certain empirical constants. In the case of heavy flavor production 
we substitute for xbj expressions for x and y from (|TT|l and (|T^ for each of the nuclei: 



and the same for the other nucleus (with x ^ y). Since Qs ~ A^^^ the saturation 
scales hke iVp^rt- In Ref. |7] dependence of Qs on A^part was calculated using the Glauber 
approach. Here we are going to use the result of that calculation and refer the reader to 
the Ref. jjj for the explicit table. 

3.2 Numerical results 

We perform numerical calculation of the charmed quark multiplicity spectrum in heavy- 
ion collisions at RHIC energies using Eqs. ()10p and ((221) • The parameters we use are: 
charm quark mass m = 1.3 GeV, the Golec-Biernat-Wiisthoff parameter A = 0.3, the 
strong coupling = 0.35 at qt < Qs and runs otherwise and Ra = 1.2 fmA^/'^ for the 
nuclear radius. In Fig. |21 we present results of the calculation at midrapidity and in the 
forward region in AA and dA collisions at RHIC energy of y/s = 200 GeV as a function of 
centrality. We observe that at midrapidity quark spectra exhibit A^'cou as expected since 



Qs{y* = 0) ~ m whereas at forward rapidities y* > 2 the Ai'part and yA^^rt scaling is 

observed for AA and dA collisions, respectively. 

We conclude that the effects of multiple rescatterings and anomalous dimension vari- 
ation do not affect the centrality dependence in the central rapidity region at RHIC 
energies for the charm quark production. However saturation significantly affects charm 
production mechanism in the forward rapidity region. 

The charmed-meson spectrum can be calculated by plugging the unintegrated gluon 
distribution ()23|) into ITUI) and convolving with the charmed quark fragmentation function 




(35) 




(36) 





d^Pix dyl 



I dz 




e z - 



HpI 



z'llJD^^{z)\ (37) 
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Figure 2: Dependence of the charmed meson yield on centrahty at midrapidity and pseu- 
dorapidity 77 = 2 for dAu and Au-Au coUisions. 
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D-rresons spectrum 




Figure 3: The charmed meson spectrum in heavy- ion colhsions at various energies. 



where g^ax = VS/2. In our calculation we use the Peterson function ^42j with e = 0.043 
inij. Since we determined that the spectrum scales with Ai'cou in the midrapidity one 
can easily infer the dA spectrum by scaling it with a corresponding number of binary 
collisions. The result is then compared to that obtained from PYTHIA event generator 
with the settings used by PHENIX collaboration jl^ (note that PYTHIA is based upon 
the coUinear factorization). In Fig. El we observe that PYTHIA spectrum is significantly 
softer than ours. Therefore, we predict much harder open charm spectrum. Of course, the 
steepness of the PYTHIA spectrum depends on the intrinsic momentum parameter ko. 
Our spectrum can probably be reproduced by PYTHIA if one takes fco — ^ ^qcd and 
fixes the ii'-factor. This large value of intrinsic would however signal the breakdown 
of the collinear factorization approach. 

3.3 Distribution in the total momentum of the heavy quark pair 

To get further insight into the dynamics of the heavy quark production it can be instructive 
to directly measure the total momentum of qq pair, which has to be equal to zero in the 
leading order pQCD calculation. Consider the relative motion of the two produced quark 
momenta in the transverse plane. In the region where collinear factorization is valid the 
transverse momentum conservation of the process gg — >■ cc tells us that the transverse 
momentum distribution must be proportional the delta function S{pi±—p2±) in the center- 
of-mass frame. Due to hadronization the delta function smears out to become a Gaussian 
of the typical width of the order of A^ independent of the energy and the impact parameter 
(centrality) of the collision. At momenta < fcgeom — Qs the collinear factorization 
breaks down and we have to use the fcy-factorization approach. The relevant process is 
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99 ~^ 99CC. Therefore momentum conservation no longer requires that the quark and 
antiquark momenta be correlated back-to-back pi± = —p2±- In general, the azimuthal 
distribution of quarks becomes asymmetric and their transverse momentum distribution is 
disbalanced. The typical difference between the quark momenta is Qs and therefore grows 
with the energy of collision and centrality. It also become larger for quarks produced far 
from midrapidity. As wee see, the direct consequence of the saturation of the nuclei wave 
functions is that the total transverse momentum of the quark-antiquark pair {pi±+p2±Y 
does not vanish. 

Generally, in a given event a number r of qq pairs is produced. Let Pi be the total 
transverse momentum of all quarks and anti quarks produced in a given event. The 
averaged total transverse momentum in events with a given r is given by 



2r \^ 2r 



d'^P±i dyl . . . d'^p±2r dy 



(3J 



2r , 



With a good accuracy quarks and antiquarks are produced in pairs, thus at this approxi- 
mation we can write 



/diV 



2 

2r 



= r{Pl),. (39) 

It is convenient to define the event-averaged total transverse momentum of the cc pair as 
follows 

2\ _ / ^7^2^ 



^Pi) = {-APDi) ■ (40) 

\ ' /all events 

We can estimate (P|) using equations (IISI),(IIlI),(II7j) together with (HH)) , (HHI) and (OT|l . 
Since in the midrapidity region of RHIC cc is produced perturbatively we find 

(Pi) ~ 0, / ~ 0, RHIC. (41) 

At larger collision energies, when Qs{y* = 0) ^ m both nuclei are saturated and we get 
(modulo logarithms) 

(Pi) ^ Ql y* ^ 0, LHC. (42) 
At large rapidities, when only one of the nuclei is in saturation we find 

(^l)^QLax- RHIC, LHC. (43) 



4 Quenching of the heavy quark spectra 
in QCD matter 

Our discussion so far has neglected the final state effects. However if hot QCD matter 
is formed in heavy-ion collisions it could strongly influence the open charm spectrum. 
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similarly to the case of light partons fOl 1201 1^ • The energy loss of heavy quarks however 
is different in one important aspect jHj: since the velocity of heavy quark is smaller 
than unity, the angular distribution of the gluons emitted in the medium vanishes in the 
forward direction ("dead cone" effect). Consequently, the coUinear singularities disappear, 
and the sensitivity of the result to the infrared cutoff is greatly diminished which allows for 
a rigorous perturbative QCD treatment even at small transverse momenta. The resulting 
energy loss is significantly smaller than for light partons. Recently this issue was also 
examined in Refs [IH] . 

The medium-modified spectrum can be written as [201 IS] 



da 



med 



da^ 



dpi 



dpi 



ip±) Qh{pa 



(44) 



provided that the energy loss is much smaller than p_L. The quantity da" /dpi_ corre- 
sponds to the quark spectrum in the vacuum. The quenching factor Qh{,P±) for heavy 
quarks has been evaluated as jUj 



Qh{pa 



exp 



L 



£iPl 
Pi. 



L 



qm 



1/3- 



+ pi, 



where 



^PA 



d 



In 



da^ 



dpi 



■{Pt) 



(45) 



(46) 



d\np^ 

q is the transport coefficient and L is the medium size. The charmed meson spectrum in 
vacuum calculated in the previous section has intrinsic momentum po — Qs- Obviously 
at p_L ^ po we have C = const and the quenching factor rapidly approaches unity. 
However, when p± ^ po we have C oc p± and the first term in the exponent in ()45|) stays 
almost constant. This implies much slower dependence of the quenching factor Qh{p±) 
on Pt- We calculate £ by substituting the quark spectrum into (HOI)- Then, using (ji^ 
we calculate the quenching factor for AA collisions (hot matter) and dA collisions (cold 
medium) which is plotted in Fig. |3 Note, that the spectrum shape in vacuum in fc^- 
factorization is much harder than that in the parton model. Therefore, the fact that 
we use /cT-factorization is both qualitatively and quantitatively essential for the result 
presented in Fig. EJ Even though Fig. E) exhibits substantial quenching for the values 
of the gluon transport coefficient that we used, it is still significantly smaller than the 
quenching factor for the light quarks computed under the same assumptions about the 
density of the medium. 



5 Summary 

In conclusion, we have calculated the charmed meson spectrum for heavy ion collisions in 
the framework of the Color Glass Condensate. We found that at midrapidity the spectrum 
at RHIC energies scales as Ai'cou, while in the forward region it shows Np3.rt{\fN^t) scaling 
in AA (dA) collisions. Our results are thus different from the predictions based on coUinear 
factorization, even after the leading twist shadowing is included. 
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Figure 4: Quenching factors for hot and cold media q = 0.2, 0.05, 0.01, 0.0025 GeV^. 
L — 5 fm. 
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Using the explicit form of the charm quark spectrum computed in this approach, we 
have calculated the quenching factors due to interaction of quarks with hot and cold 
media. Although at very high the quenching factor approaches unity, in a wide region 
of 2 < < 15 GeV it stays almost flat. This feature of the quenching factor is a 
direct consequence of the saturation of the nuclear wave functions which brings in the 
dimensionful scale po ~ Qs- 

We have suggested to determine the saturation scale by measuring the total transverse 
momentum of charmed mesons produced in a given rapidity interval. We expect that 
unlike in naive leading order perturbation theory, the total transverse momentum of the 
pair will not vanish. On the contrary, it will grow with the energy, rapidity and centrality 
in the same way as the saturation scale. 
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Appendix 

In the Appendix we cite the heavy quark production amplitude A{gg ggcc) as given 



m 



-^{gg ggcc) 



2iV, 



2{N^ - 1) 



A 



nab 



where 



A 



ab 



t — m? u — m? 



A 



nab 



2 

+ 



1 



[t — m?'){u — m?) 



S sxy 



+ 



S ^2 X2yi 
2 t-m^ 



s \t — 
A 



\x y 
A- 



A = S 



Sxiy2 - Sx2yi + qj^— - q^^— + -(w-t) + ^iql± - ql±) 



u — 
1 



x y 2 2 



^ = (?1 - P2) , u = [qi-pi) , S = M = {pi + P2 

qi = ypA + qi±, q2 = xps + g2± 

Pl,2 = yi,2PA + Xi,2Pb + Pl,2±, 

where 2pA ■ Pb = S. 



(47) 



(4J 



(49) 

(50) 

(51) 
(52) 
(53) 
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